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Abstract: Using the AdS/CFT, we establish a correspondence between the 
intricate thermal phases of R-charged AdS^ blackholes and the R-charge sector of 
the N=4 gauge theory, in the large N limit. Integrating out all fields in the gauge 
theory except the thermal Polyakov line, leads to an effective unitary matrix model. 
In the canonical ensemble, a logarithmic term is generated in the non-zero charge 
sector of the matrix model. This term is important to discuss various supergravity 
properties like i) the non-existence of thermal AdS as a solution, ii) the existence 
of a point of cusp catastrophe in the phase diagram and in) the matching of saddle 
points and the critical exponents of supergravity and those of the effective matrix 
model. 
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1. Introduction 

The AdS/CFT correspondence |]I[] implies that the phases of string theory can be 
studied by studying those of the dual gauge theory ([0, [5, Q). In the case of 
type IIB string theory in AdS^ x the phases can be studied using the large 
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N limit of a unitary matrix model(P, The unitary matrix is the finite 

temperature Polyakov loop which does not depend on the points of S^. This fortunate 
circumstance is due to the fact that in the Hamiltonian formulation, N=4 SYM theory 
at a given time slice, is defined on the compact space and the SO (6) scalars are 
massive because of their coupling to the curvature of 5*^. These facts imply that, in 
principal, one can integrate out almost all the fields and obtain an effective theory 
of the zero mode of the gauge potential Aq. 

Using this method a detailed correspondence of the critical points of the gauge 
theory effective lagrangian and the critical points of supergravity (discussed by 
Hawking and Page ||1 1|| ) can be constructed at the leading order of the expansion. 



These are AdS^ and the small and big black holes (as [|T^ we refer to these as SSB 
and BBH.) It turns out that in the gauge theory these critical points are in the 
gaped phase, where the density of eigenvalues vanishes in a finite arc of the circle 
around which the eigenvalues are distributed. The closing of the gap, corresponds 
to the Gross-Witten(GW) phase transition([|l^, |l5l). In a window around this 
transition, the supergravity description of string theory is likely to smoothly cross 
over into a description in terms of heavy string modes. ^ 

In this paper we extend the discussion of the correspondence between R-charged 
AdS^ blackholes(|]T9|, [2^, |18|, 0), and the effective unitary matrix model(|l3|, 
T^, [T^, 0). R-charged black holes are known to have a rich phase structure in the 



canonical and grand canonical ensemble. In the canonical ensemble the fixed charge 
constraint, contributes an additional logarithmic term log(Tr?7Trf/^) involving the 
order parameter, to the gauge theory effective action. This term is crucial for 
matching with supergravity. We analyze the implications of this term in the large 
N limit and compare with the various supergravity properties like the existence of 
only blackhole solutions in the canonical ensemble and also the existence of a point 
of cusp-catastrophe in the phase diagram. 

The plan of this paper is as follows. In section |^ we give a brief review of 
charged AdS^ blackholes. In section |^ and section ^ we discuss the effective action 
of the gauge theory at zero and small coupling, in the fixed charged sector. At zero 
coupling there is exactly one saddle point and the value of TrUTrU~^ at the saddle 
point is always non-zero. For a small positive coupling there are two stable and one 
unstable saddle points, all with a non-zero value of TrUTrU^^ . They merge at the 
GW point. In section (|^) we discuss the model effective action at strong coupling. 
Here too, there are three saddle points,two stable(I,III) and one unstable(II). In the 
region p > |, / and /// can be identified with a stable small blackhole and stable 
big blackhole respectively. Saddle point // is identified with the small unstable black 



^Phases of large N gauge theory is also discussed in |l^ 

^In ||l2|, besides the vicinity of the Gross- Wittcn transition, authors also studied corrections 
and presented formulas for the partition function in the vicinity of blackhole nucleation and the 
Hagedorn transition. 
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hole. The merging of saddle points leads to critical phenomenon whose exponents 
can be calculated and shown to agree with supergravity. This is discussed in section 
||. We have also calculated the o(l) part of the partition function near the critical 
point. 



2. R-charged blackholes in ^4^55 and critical phenomena 



The R-charged AdS^ black hole and relevant phase structure were discussed by 
A. Chamblin, R. Emparan, C. V. Johnson and R. C. Myers ( |jl8[ , ||T^ ) . Here we 



review their result. The Einstein-Maxwell-anti-deSitter (EMadS^+i) action may be 
written as 



d 



-9 



R-F' + 



n[n 



1) 



i?2 



(2.1) 



where R is the Ricci scalar and R is the characteristic length scale of AdS. The 
metric of the Reissner-Nordstrom-anti-deSitter (RNadS) solution is given in static 
coordinates by 



ds^ = -{l-^ + ^ + ^)de + "^^^ - + r^d^^ 1, (2.2) 

The parameter q is proportional to the charge 

Q=v/2(n-l)(n-2) (^) g (2.3) 
and m is proportional to the ADM mass M of the blackhole. 

M = ("-^^""-^ m (2.4) 
ujn-i is the volume of the unit (n— l)-sphere, and the gauge potential is given by 

^o=\-^j^;h + ^\ (2-5) 

where $ is the electrostatic potential difference between the black hole horizon and 
infinity. 

For n=4 the solution {^^T) can be considered as a rotating black hole in AdS^ x 



with angular momentum in the internal space S^.^ The symmetry group of 5*^ is 



^The EMadS system described here may be thought as the dimensional reduction of AdS^ x 
to AdS^. Generally introducing angular momentum in will distort S^. This distortion is not 
taken into in the EMadS reduction j22| We thank S.Trivedi for pointing this out to us. However 
these details do not change our main consideration. 
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5*0(6) and the black hole we are discussing has equal U{1) charges for all the three 
commuting U{1) subgroups of S0{6), the R-symmetry group of the N = 4 SYM. 
Hence we are dealing with a system which has the same chemical potential /i for all 
three U{1) charges in the grand canonical ensemble or equivalently three fixed equal 
U(l) charges in the canonical ensemble. 

2.1 Equation of state 

In order to discuss the thermodynamics, we consider the Euclidean continuation 
(t^ir) of the solution, and identify the imaginary time period f3 with the inverse 
temperature. Using the formula for the period, /5=47r/l^'(r+) (for a review see p^), 
we get 



This may be rewritten in terms of the potential as: 



(n 



2)Z2(l-c2$2) + nr2 



[2.6) 



[2.7) 



The condition for euclidean regularity used to derive ( p.6|) is equivalent to the 
condition that the black hole is in thermodynamical equilibrium. The equation (|2.6| ) 
may therefore be written as an equation of state T=T{^,Q). From this equation 
of state we see that for fixed $ we get two branches, one for each sign, when the 
discriminant under the square root is positive |]18[. For fixed Q, T($) has three 
branches for Q<Qcrit (Let us call them /, //, ///) and one for Q>Qcrit- The critical 
charge is determined at the "point of inflection" by {dQ / d^)j, = {d'^Q/d^'^)rp=0. 

The qualitative features of /5(r+) for varying q are shown in Fig |l]. There is a 



l/T 



Figure 1: Plot of P{r+) for q increasing from the left. The third graph from the left is 
for Qcrit 



critical charge, gcrit, below which there are three solutions for r+ for a range of values 



of T, corresponding to small (/), unstable(//), and large (///) black holes. For fixed 
Q < Qcrit only branch I is available at low temperatures. At T = To2(g), there is a 
nucleation of two new solutions, II the unstable small black hole solution, and III 
the stable big black hole solution. As the temperature is increased further, the black 
hole // approaches black hole / and at T = To2{q) the two solutions merge. 

As q is increased further, Tq2 increases, whereas Tqi decreases. At g = qcru we 
have Tcrit = Toi = 7o2- At g^rii and Tcru all three solutions merge. In the language 
of catastrophe theory this is a cusp catastrophe. As we increase q beyond qcrit there 
will be just one solution for all temperatures T. 

We wish to take note of some properties of the phase diagram. 

1. Thermal AdS is not a solution and all three branches of the solution represent 
black holes. 

2. There exists a critical point where the three solutions of the system merge. It 
is a point of fold catastrophe. 

2.2 Critical Phenomena 

The critical point {Qcit,Tcrit) may be approached from various directions in the 
parameter space. If we set T = Tcru, then the equation determining (r — rent) 
takes the form {rcru is the radius of the critical black hole) 

(r - r,„t? = C{Q - Qcrit) (2.8) 
C is a numerical constant. The critical exponent here is ^, since 

{r - Vcrit) 0C{Q - Qcrit)^ (2.9) 

As discussed in [|18| (Figl6), the critical point may also be approached through 
the coexistence line in the parameter space. The coexistence line is the line with the 
property 

Si = Sin (2.10) 

for the parametric range q < qcrit- It is the line where the Hawking-Page (first order) 
transition from the small black hole to the big black hole takes place. As we approach 
the critical point through this line, we have the relation 

(r,-r,,) oc (T-T,„t) 5 (2.11) 

In the following we will present an understanding of these properties. Before we 
do the matching with supergravity we would like to present a discussion of the gauge 
theory in the limit when X = g'^N = and also when A << 1. In these cases we of 
course can not compare with supergravity which requires A >> 1. 
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3. Free YM theory 



3.1 Effective action with chemical potential 

In this section we briefly review (see 3,0) the effective action for a free SU{N) 
Yang-Mills theory (with adjoint matter) on a compact manifold in the large N 
limit. The basic idea is to integrate out all fields in the theory except for the zero 
mode of the Polyakov line. The partition function is then reduced to a single unitary 
matrix integral. 

Expanding all fields in the gauge theory in terms of harmonics on 53, the theory 
reduces to a zero dimensional problem of free N x N Hermitian matrices 

C = lj:Tr[{D,M^r-u;lM^] (3.1) 

The sum in ( ^.l] ) )is over all field types and their Kaluza-Klein modes on S3. u;„ is 
the frequency of each mode. The covariant derivative in ( ^.ip is 

DtMn = dtMr,-t[Ao,M^] (3.2) 

Aq comes from the zero mode (i.e. the mode independent of coordinates on S3) of 
the time component of the gauge field and is not dynamical. The partition function 
of the theory at finite temperature can be written as a unitary matrix integral by 
integrating out all fields in ( |3.1| ) except for Aq. Hence we have 

Z = JDU n„ (det^rf, (l - e„e-^'^"f/)) (3.3) 

where U = exp{ipAo) is a U{N). detadj denotes the determinant in the adjoint 
representation and e„ = 1 (—1) for bosonic (fermionic) M„. The above equation can 
be expressed as 

Z = f DUexpiY ^^)TrU^TrU-'^) (3.4) 

where 

zn{[5) = ZB{nP) + {~ir+'zF{nP). (3.5) 

Here zb{i3), zf{i3) are the single particle partition functions of the bosonic and 
fermionic sectors respectively ( see ^ for the explicit formulas for z{f3) in various 
gauge theories). 

If we introduce a chemical potential n = (logm), the formula for the partition 
function changes to 

ZIP, /X] = / Z}f/e(S:^i ^TrU-TrU--) ^3^g^ 
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where 

Zn{(3, fi) = ZB{nf3, nfi) + {-l)''^^ZF{n(3, nfi) (3.7) 

zb{P, fi), zf{P, ^) are now the single particle partition functions with chemical 
potential fi. Hence 

ZB{f3,fi) = ''M-EiP + Qifi) (3.8) 

bosons 

ZF{(3,fi)= exp{-Eif3 + Qifi) (3.9) 

fermions 

Qi is the charge of the state whose energy is Ei (see 0). 

As we have already mentioned, we wish to describe a system which has the same 
chemical potential fi for all three U{1) charges of the R-symmetry group S0{6) in 
grand canonical ensemble. Equivalently, we can work with fixed and equal values of 
the U{1) charges in the canonical ensemble. 

Let us for simplicity confine ourselves to the bosonic sector of the = 4 SYM 
theory. The gauge fields have no R-charge. The six scalars (pi {i = 1, ...,6) can be 
grouped in pairs of two. We define 

0+ = 0i + 0- = 0i - ^3^^Q^ 

(We can similarly define complex fields for the other two pairs.) 0^ have charge ±1 
for each of the three commuting U{l)s of 5*0(6). Hence, if we consider the single 
particle partition function for these fields, it will be 

z[x, /i] = {exp{+fi) + exp{—fi))z[x, 0]/2 = cosh(/i)z[a;, 0] (3-11) 

where z[x, 0] is the single particle partition function without any chemical potential. 

3.2 Canonical Ensemble 

We will now discuss the free gauge theory partition function for a canonical ensemble 
with constant charge Qq, by introducing a delta function S{Q — Qo) in the functional 
integration of the gauge theory. Q = Q[(j)] is the corresponding functional for the 
charge which we want to keep fixed. In our case Q is just the functional for R-charge 
in gauge theory. 

The fixed charge partition function is defined by 

ZW,Qo) = J DXelo^ms^Q _ g^) (3.12) 

= J DXe/o"^!^] J exp{tfiQ)exp{-ifiQo)dfx (3.13) 

= J dfiexp{-ifiQo){J DXefo^^^^e'^"^) (3.14) 

= J dnexp{-ifiQo) J DUexp{i:Zn[l3,iij]TrU''TrW') (3.15) 
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where z„[/3, z/i] = zl[(3, 0] + cos(/i)zf 0] + cos(f )zf 0]. 

We can now make the approximation that "i/i] | for n > 1 is neghgible 

in comparison to Neglecting the contribution from the n > 1 modes we 

arrive at a model which contains only TrUTrU~^. Using the specific formula for 
z[x,fi], of the bosonic ^ sector , 

Z{/3,Qo) = J (iyuexp(— z/iQo) J DUexp{{a + ccos{fi))TrUTrU~^) 

= J DUexp{aTrUTrU~^) J (i/iexp(— i/i(5o)exp(ccos(/i)Trf/Trf/~^) 

= J DUexp{aTrUTrU-^)lQ,{cTrUTrU-^) (3.16) 

Here a{P) = z^[P, 0],c(/3) = z^[P, 0] and for convenience we did not show the explicit 
P dependence in the equations. /„(x) is the Bessel function. 

Hence we end up with a matrix model with an effective potential 

Seff = a{TrUTrW) + \og[lQ,icTrUTrU-')] (3.17) 

where a > 0, c> 0. We define = TrUTrlJ-'^/N^ to get 

S^ffip) = N\ap' + il/N') log[/Q„(ciVV^)]) (3.18) 

It may seem that the logarithmic term is suppressed by a factor of l/N"^ and 
hence negligible for large N. But this need not be the case because, in the semicalssical 
large N limit, we must deal with a system with a charge of order N"^. Hence we define 
Qq = N^q {q ~ o(l)). Using the asymptotic expansion of Inijix) for large n, the 
effective action ^ becomes 



SefAp) = N\ap' + qd{l + ^y) + log( 1^=^))) + 0(1) (3.19) 

3.3 Phase Structure 

To understand the phase diagram of this model at large N, we have to locate the 
saddle points of ( |3.19D after including the relevant contribution from the path integral 
measure depending on whether p < ^ or p > ^7 

^This approximation can be thought of as a low temperature approximation. This is because 
at low temperatures z^[f],0] approaches zero as e~'^". Hence the higher z„ are suppressed relative 
to zi. It is also true that for all temperatures, z„ < zi and for very high temperatures we have 
Zn ^ As an example, the total contribution for all other z„, even near hagedorn transition in 
free N=4 SYM theory, is only about 7% of 

Unitary matrix models involving TrC/" , n > 1 has been discussed in ||2^, . 

^Effect of the fermions is discussed in appendix 

^It should be noted that when Q=0 wc should use the asymptotic expansion of Io{x) for large 
X. Then we get the expected answer Seff — (a + c)TrUTrU^^ which is same as a model without 
any constraint on charge. 

^The term in the right hand side of equation ( |3.21 ) originates from the path integral measure 
over an unitary matrix. ( see p7[ , Appendix of |l3|). 
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Differentiating Seff{p) we get 



a + b 



.2 



a + 4(l + V)^+0(VQ) (3.20) 



Hence the equations to solve are 

ap + -{l + —p )2 = p,p<- 
p q ^ 



+ + = __,,>_ (3.21) 



The left side in ( p.21|) can be written as 



ap + ^(1 + 4p')^ = ap + cp + ^- — (3.22) 

So fixing the charge gives rise to a term of type ^-^-t which has some 

important properties. 

1. For all values of g > 0,c > 0, this term is a decreasing positive function of p 
and it diverges as p — > 0. 

2. For all values of c > it is a monotonically increasing function of q. 

We can now discuss the solution of this model at = oo. Let us assume that 
we are discussing the phase where a{T) + c(T) < 1. This condition is valid for low 
temperatures since a(T), c(T) ^ as T ^ 0. It should also be recalled that without 
any charge fixing the hagedorn transition occurs when a{T) + c(T) = 1. Unlike the 
situation with no charge, here we have a function, on the left hand side of (|3.21| ), 
which diverges as p ^ 0. Hence p = can not be a solution. We get only one 
solution at a finite value of p which we will describe in the next paragraph. 

Equation ( |3.21| ) is solved in the region p < \ with solution 



4 

^ " (1 - a)2 - c2 " (l-a-c)(l-a + c) ^^'^^^ 
The self consistency condition for a solution in the region p < | is 

p' = n — T2 2 < 4 (3-24) 

1 — ar — c"^ lo 
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At low temperatures the condition is naturally satisfied for a small enough value 
of q. If we gradually increase the temperature (i.e. the value of a(T) and c(T)) 
while keeping the value of the q fixed, then the value of p at this saddle point will 
increase. At some temperature T^iq), p will become equal to ^. Since the measure 
part (i.e. right hand side of (|3.21| ) ) has a third order discontinuity at p = |, we will 
get a third order phase transition at the temperature T3. From ( |3.24D we have the 
following condition at T3 



q^ 1 



ar — (? 16 



(3.25) 



If the temperature is increased beyond T3 then we have to solve ( p.21|) in the region 



p>\. 

If we increase g, then the value of p at the saddle point for a fixed temperature 
will increase. At some ga we get a third order phase transition satisfying 

16g3(T)2 = (a - 1)^ - (3.26) 

Since the minimum value of a(T) and c(T) is zero, the maximum possible value of 
g3(T) is g^^jj — If we increase the g beyond qcrit-, the saddle point will always be 
confined in the parameter range p> \- Consequently as we increase the temperature 



1 

4' 



from zero we will not get a third order phase transition for g > qcrit 

This free model, at zero gauge coupling {Ig » R in bulk), has some similarities 

with AdS^ black holes in a fixed charge ensemble. However unlike the three black hole 

branches in AdS^, we get only one branch in the free theory. But most importantly 

the solution always has a nonzero value of p. 

It should be recalled that before the Hagedorn transition, a free gauge theory 

with zero charge has the solution p = [|, 0]. 

Some properties of the free theory will be important when we analyze the 

situation for the weakly coupled gauge theory. Just above the temperature T^iq), the 

difference of the two sides of ( p.21| ) can be expanded in the region p > |. Defining 

p = i + x,(x > 0) the difference is 

-e(g)x - Cix^ (3.27) 

Here e(g) > and e(g) — > as g — 0. It is important to note that Ci > 2 because 
the measure function (i.e. right hand side of ( |3.21|) has a third order discontinuity 



at the point p = \- We will also discuss the significance of this in what follows. 
4. Small coupling model 

We will now discuss the problem with a small non-zero gauge coupling A = Qym^ ■ 
By AdS / CFT correspondence it corresponds to a finite string length in AdS. It has 
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been shown in 



that by considering a phenomenological model of type 



S[TrUTrU-^] = a(A, T){TrUTrU-^) + ^^^^P {TrUTrU-^f , b>0 (4.1) 

we can map out the possible phase diagram of type II B string theory in AdS^.^ Even 
though the model in ( |4.1| ) can be derived from a weak coupling analysis of the gauge 
theory, it can be thought of as a phenomenological model describing supergravity in 
AdS^. 

We are motivated to discuss the fixed charge ensemble in the same spirit. Let 
us add a small interaction term to (|3.18|) . The effective action is then given by 



c J2 



S^jfip) = N\ap' + bp' + qJ (1 + + log( f=^))) + 0(1) (4.2) 

V 1 1 + y'l + 

Here b is proportional to A and is also a function of charge. Depending on the theory 
considered, the sign of b can be either positive or negative. It has been shown in 
that b is positive for a pure YM theory. In the following discussion we will assume 
that this is the case in order to motivate a similarity with the supergravity picture. 
The equations determining the saddle points of (|4.2| ) , including the contribution 
from the path integral measure, are 



(a + c)p2 + 26p* + 
(a + c)p2 + 26p^ + 



4 , 1 ^2 ^ 1 



P 



+ + 4(1 -p) 



(4.3) 



In what follows it will be useful to introduce a function 



M(p) = p^p<- 



4(l-p)'^>2 



(4.4) 



M(p) is an increasing convex function and is the right hand side of (O). It is also 
useful to introduce D{x) = S'^jj-{x)~M{x). Eqn ( ^Tsj) is then equivalent to D{p) = 0. 



It has been shown in that for an interacting model with zero charge, 

and 



l-2w 



we get nucleation of black holes along the curve given by a = (i_^,)2(i^^,) 
^ ~ {i-w)^i+w)s , 1 > w > 0. Here we want to analyze a similar type of phenomenon. 
Let us consider the different cases. 



^In fact in [n2[ an arbitrary convex function is considered, and shown to map out the phase 



diagram of IIB theory in AdS^. The simphfied model (4.1) leads to similar qualitative result 
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T is varying and q is small 



Let us start with a value of charge which is small (i.e. g <C 1) and let us increase 
the temperature from zero. At low temperatures all the parameters a, c will be small 
(for small T these parameters have a dependence like e~^^, c is a constant ). Hence 
we get just one solution for p < | which we call /. There is no solution for p > | 
(see the topmost curve of Fig ^) because the left hand side of ( |4.3| ) is less than the 
right hand side. 

The situation is quite similar to supergravity where for small charge and low 
enough temperatures we get a stable small black hole solution. ^ 



0.001 



-0.002 




Figure 2: Plots of D[p) with increasing a = c from the above and with fixed h and q < qcrit 

The function M(p) (i.e. right hand side of ( |4.3| )) is a convex increasing function. 
Hence we will generate new solutions of (^4.3| ) in the region p > | as we increase 
temperature (i.e. a(T), c(T) as discussed in appendix ^) keeping q fixed (Fig Fig 
^). The new solutions will always come in pairs (Fig |^). Let us call the solution 
nucleation temperature as Toi(g). At T = Tqi we have 

D{p) = S',jfip)-M{p)=0 

D\p) = (4.5) 

As we increase the temperature further (i.e. ci(T), c(T)) the function on the left 
hand side of ( [4.3|) will also increase. Hence the solutions will start to separate. Let 
us call them // and ///. Here pni > pii- Also, /// is a stable saddle point whereas 
II is an unstable one. They are similar to stable big and unstable small black holes 
in supergravity. As the temperature is increased beyond Tqi the value pn decreases 
whereas pm increases. At some temperature T^iq), we will have S{III) < S{I) 
and consequently we expect a first order phase transition. At a temperature T^p, 

^We should keep in mind that p < ^ is not the supergravity regime. The solutions of gauge 
theory effective action there should be represented as excited string states. 
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Figure 3: Top three graphs of Fig 2 showing emergence of two saddle points in the region 
P>k 

the dominant saddle point of the system changes from I to III. As the temperature 
increases the saddle point // goes through a third order transition when p{II) crosses 
the P = \ point. Call this temperature T3 which is determined by the following 
relation between the parameters 

(a+^)+4g(l + ^)^ = l (4.6) 

Increasing the temperature further, the saddle point // approaches the saddle 
point I and they merge at a temperature T02 (4t/i graph from above in Fig ^ ). In 
the language of catastrophe theory this is a fold catastrophe. For T > T02, the only 
saddle point is ///. This then is the thermal history as we increase the temperature 
for a small q. 

In summary at low temperatures, we have only one saddle point / and then two 
new saddle point //, /// are created at Tqi. As the temperature increases the saddle 
points /, // merge at a temperature T02. Beyond that we have only one saddle point 
///. In the next paragraph we will discuss what happens when we increase the value 
of the q. 

Varying q: 

Let us discuss how the various temperatures discussed above change as we increase 
the value of q. 

1. The first is Toi(g), the nucleation temperature for saddle points // and ///. 
Tqi will decrease as we increase the value of q. This is so because all three 
terms in the right hand side of (|4.3|) are positive and increasing functions of p 
and the left hand side is a positive convex function. Consequently the saddle 
point value of p at Tqi will also decrease. 
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2. The temperature T02 at which the stable saddle point / and unstable saddle 
point II merge and the value of p at T02, will increase as we increase q. The 
reason is that all the coefficients in ( [4 .31) are positive and hence increasing the 
value of q increases the function on the left hand side. 

As we increase the value of q further, Tqi and T02 will become equal for some 
value of g = gcrit- In the language of catastrophe theory this is a cusp catastrophe. 
Corresponding to qcrit there will be a Tcrit = Tqi = T02 and a value of pcrit (saddle 
point p at qcrit , Tcrit)- As we increase the q beyond qcrit, we do not get any new 
saddle point and consequently there is only one saddle point for all temperatures. 
We will discuss the physics near this phase transition in detail in what follows. 

Value of Pcrit • 

Let us determine the value of pcrit- As already discussed at the end of the previous 
chapter, in ( ^.27| ) C is always a finite quantity. Hence, ( |4.3| ) cannot have three 
solutions in the region p < | for small b. Therefore the saddle point / is always in 
the region p < |. Whereas the saddle point /// will be in the region p > |. Hence 
the only place where these three saddle points can meet is pcrit = | which is also the 
point of the third order phase transition. 

Physical Interpretation 

Before proceeding further we will briefly discuss the bulk interpretation of saddle 
points of weakly coupled gauge theory. Weak coupling in gauge theory means 
Is ^ RAdSa- Hence the supergravity picture is not valid in the bulk. However at 
large N, the string coupling (i.e. -^) will be small and we may conclude that the 
saddle points discussed above can be described by exact (in all orders in Ig) conformal 
field theories. These CFTs are characterized by the values of q and p at the saddle 
points. 

We would like to end this section by emphasizing that the coincidence of the 
three saddle points at p = |, is a property of the weak coupling (A << 1) hmit of 
the gauge theory. In what follows we will see that this fact is not necessarily true 
at strong coupling. We will show that there the coincidence happens in the gaped 
phase where p > |. 

5. Effective action and phase diagram at strong coupling 

In this section we will discuss the effective action and the phase diagram in the 
strongly coupled gauge theory which is dual to the supergravity (discussed in section 
1^) regime of IIB string theory. 
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5.1 Finite temperature effective actions in the gauge theory 

Let us first summarize the situation in the zero charge sector. Tlie propagator of 
adjoint fields in the free gauge theory, on a compact manifold S^, coupled to a space 
independent Aq is given by (see ||29[| ) 

n=oo 

G'^,,ix,t,y,0)= 5: (f/")}Go(x,t + n/3,y,0)(t/-")J (5.1) 

n=—OD 

where Go{x,t,y,0) is the zero temperature Green's function and U is the constant 
Polykov line. We know that at any temperature Go{x, t + n(3, y, 0) > Go{x, t, y, 0) 
and also at low temperatures 

Go{y,t + n(3,x,0) ^ e-""^ (5.2) 

Using the above Green's function one can develop the large N diagrammatic to 
arrive at an effective action involving Zn invariant terms built out of products of 
trU^ . In fact one can imagine integrating out all the modes trf/" for n > 1 in favor 
of trUtrU^^. In this way one gets an effective action of the form 

n=oo TrUTrU 

Seff= E «n(AA)( ^ ^ )" (5.3) 
n=l 

As one increases the coupling constant we would expect that the form of the 
effective action would remain the same except that the dependence of the parameters 
on temperature and the 'thooft coupling would change. 

5.2 Non-zero charge sector 

If we include the fixed charge constraint, as in ( |3.15| ), then we get the following 
expression for the fixed charge path integral 

Zg = J dfie'"^ J Df/e^'(S"=o°°^"('''^'^)=°^("^)) (5.4) 

For large N we can do the n integral by the saddle point method. The saddle point 
of n is on the imaginary axis. Hence we set fi = im, to get the saddle point equation 

n=oo 

q= nSnjp, A, /3)sinh(nm)), q = Q/N'^ (5.5) 

n=l 

At small values of p, S'„(p) goes as p" and hence in the p ^ limit we can 
approximate the equation (|5.5| ) as 

q ~ Cpsinh(m) (5.8) 



"'^^ Temperature here is nieasured in units of -it—. 

^^Introduction of a chemical potential changes the formula ( |5.l[ ) as 



n—oo 



G^^%{x,t,y,0) = e'^ J2 {Un]G^{^,t + n(3,y,0){U-^^)'; (5.6) 



n— — OO 
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where C is a constant independent of p. The solution is 

m ^ glog(p) + C,p ^ (5.9) 



Substituting m in ( |5.4| ) we get a logarithmic term for p. We conclude that the 
logarithmic term is a general feature and it implies among other things that TrU = 
is never a solution in the non-zero charge sector. 

5.3 Model effective action at strong coupling 

Following our previous discussion, we will include the generic logarithmic term in 
the effective potential for a fixed non-zero charge, and proceed to analyze the saddle 



point structure following [112 



Our proposal for the gauge theory effective action is 

5, = 5(a(T),6(T),...,p) + glog(p) (5.10) 
and the saddle point equations are 

pFiaiT),biT),ciT),...p') + q = p\p<^ 
pF{a{T)MT). ciT), ...p') + q = P > ^ (5.11) 

Where F{p) = S'^jj-{p). We assume that 

1. F{x,T) is a monotonically increasing function of x and F(0,T)=0 

2. Value of F{x,T) increases for fixed x as we increase the temperature and 
F(x,0)=0. 

These global properties of F{x) reproduce the phase diagram of supergravity. 
Analysis of solution structure 

Let us consider the function D(T,p) = pF{p,T) — M{p) (Where M(p) is the 
contribution from measure appearing at the right hand side of ( |5.11|) ). At T = 0, 
F{p,T) is zero. Hence D(T,p) is a monotonically decreasing function of p at T = 0. 

We know that at T = Tqi a pair of two new saddle points appear at poi > |. 
Hence at T = Tqi we have D(Toi,poi) = and D'(Toi,poi) = 0. D(Tqi,p) has a 
zero for p = and it is a decreasing function in the neighborhood of p = 0. It again 
increases and become zero at poi and then the function again decreases as p — > 1. 
This implies that the function has a local maximum and local minimum. 

In summary: 

where 

G^(x,t + n/3,?;,0) =e'"^Go(x,< + n/3,y,0) (5.7) 

Hence in each order in perturbation theory the terms containing cos(n/i) also get multiphed by 
{TrUTrU^^)™', m > n or the higher operators like TrU" which can be integrated out to give again 
a term like {TrUTrU-^)"" . 
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D 




Figure 4: Plots of D{T, p) for T = 0, T = T^riuT = Tqi from below. 

1. D{0,p) is a monotonically decreasing function of p. 

2. D{Tqi,p) has a maximum and minimum. 

There is a temperature T^ru at which the local maximum and local minimum 
appear (Fig|5.3|). Let us call this temperature Tcru- At Tcru the curve D(Tcrit,p) will 
have a point of inflection at p = pcrit, say. Let the value of D{Tcrit, Pcrit) = Qcrit- 

Increasing the value of q from zero we need to solve the equation D{T,p) = q. 
We will get a solution for a non-zero value of p. Denote this solution as /. As the 
temperature increases, two new solutions appear at T = Tqi. Call the stable solution 
as ///, and the unstable solution II. As the temperature is further increased to T02, 
the unstable solution // and the stable solution / merge. For T > T02 , the only 
solution is ///. 

As q approaches qcrit from below the two temperatures Tqi and T02 approach each 
other. At q = qcrit, we have Tqi = T02 = Tcru- If we increase q beyond qcrit only one 
solution appears for all temperature. These facts are consistent with supergravity 
solutions (section 

With a sufficiently sharp rising function F{T, p) in ( |5.11| ) we can obtain this 



critical point in the region p > ^. . As the function D{p,T) is smooth in the region 
p > I the second derivative of the function D{T,p) will vanish at the inflection 
point, and we will get a third order phase transition. We can calculate the partition 
function in suitable double scaling limit near the critical point. This is discussed in 
section ||. 

A specific example 

We will now illustrate the above phenomenon in a simple model defined by 

F{p) = ap + bp^ (5.12) 



^^Which describes supergravity in the bulklll 
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where a,b > 0. 

We will determine the parameter ranges of a, 6 for which all the three saddle 
points of ( |5.12| ) are in the range p > ^- 




91(.PF{P) - JTJ^^ > (5.13) 

Putting the value of p = | in the above inequahty we get the following constraints 
on the parameters a + b < 1 and a + 36 > 2. Simplifying we have 6 > | and a < ^. 

As the coupling become stronger, we expect that b is not necessarily small and 
will be of o(l) or greater. All the saddle points of ( |5.11| ) are then naturally shifted to 



the region p > 2- Here, as was discussed in [|T2|, we can expect to match the solutions 
of the gauge theory with those of supergravity. The stable saddle point I corresponds 
to the stable black hole branch / of supergravity. And unstable saddle point // is 
matched with the unstable black hole branch II. The stable saddle point /// is 
matched with the big stable black hole in supergravity. With this identification 
the thermal history and critical behavior of the gauge theory, discussed earlier in 
this chapter, match with the thermal history and critical behavior of supergravity 
(discussed in section and ^ 



6. Universal neighborhood of critical point and the critical 
exponents 

Let us consider the effective action Stot{p,T,q) which includes the contribution 
from the path integral measure over an unitary matrix. The derivative of Stot 
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with respect to p, say G{p,T,q), gives the saddle point equations ( |5.11| ). We have 



already discussed that by a suitable choice of parameters the critical point appears 
in the region p > |. This critical point is a third order critical point because three 
saddle points of the system merges here. Hence the first and second derivatives of 
Gtot{p,T,q) with respect to p vanish at p = Pcrit,q = QchuT = T^ru- Expanding 
G{p, T, q) around the critical point , we get 

G{p,T,q) = {Spf^ + {ST)dTG+iSq)d,G+{Sp)iST)dpdTG+iSq){ST)d,d,G (6.1) 

Let us fix T = Tcru or 6T = 0. Then the equation ( |6.1| ) has one solution. In 
order to know how the saddle point value of p approaches pcrit (^P 0) as 6q — > 0, 
we equate the leading part of ( |6.1|) to zero. 

Q3q 

{5pf^ + i5q)d,G = (6.2) 



Hence dp (x 6q'i and we get the same universal exponent \ as in supergravity [ff8 



6.1 Partition function near the critical point 

Near the critical point we can write the Stot as 

Stot = StotiPcrit, T,ru, qcrit) + (^p)'^ + {^0)9,8 + {5q){6 p)d ^8,8 + 0(5p5) (6.3) 

1 3 

If we define a double scaling limit p = x, N^q = 2; , we can write the o(l) part of 
the partition function, after suitable rescaling of the variables, as 

Z2 oc / c/xe-(^'-"^) (6.4) 



This can be calculated in a power series 

6.2 Approaching the critical point through a line of first order transitions 

Another type of double scaling limit is possible in this problem. We can set^^ 

{5T)dTG + {5q)dqG = (6.6) 



by choosing a suitable relation between 6T and 6q. Using ( |6.6| ) in ( |6.1| ) we get 



(5p)'^ + {6p)i{6T)d,dTG + {5q)d,d,G) = (6.7) 
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It is same as following the HP (first order) transition line in parameter space. 
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with the solutions 



6p = 



(6.8) 



We can expand Sejf as 



Seff ~ Scrit + (M 




+ Ci{6T){6py + OT 



(6.9) 



4! 



where OT are terms independent of 6 p. Defining a suitable double scaling limit. 

3 1 

N2ST = z, NiSp = X and a suitable rescaling of the parameters we can evaluate the 
0(1) factors in the partition function as 



where z < 0. 

7. Conclusion 

In this paper we have studied the logarithmic matrix model generated by fixing the 
R-charge in the gauge theory partition function. In the free gauge theory it has been 
shown that there is no solution with p = {AdS type solution). We then studied the 
effect of adding an interaction term in our model and discussed the generic nature 
of the logarithmic term even at arbitrary value of the coupling. We identified the 
supergravity saddle points and their critical behavior which was discussed in ([0). 

Our main aim was to give another example of the utility of unitary matrix 
methods in providing a non-perturbative dual description of blakholes in AdS and 
to understand the relation between matrix models and string theory in general. It 
would be interesting to consider an effective unitary matrix model to describe phases 
of Kerr-Ads black holes. 
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oc / dxe-(^'+2..^) oc £ ^r(^ + ^) oc ^ze'^K.i^) (6.10) 



n=0 
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A. Appendix: Inclusion of Fermions 



Including the contributions from the fermions of = 4 SYM theory change ( p.l6| ) 
to 

Z{p,Qo) = J DU J dfiexp{N^{a + ccos{fx) + dcos{^))p'^ - ifiQo) (A.l) 

Where d{f3) is the single particle partition function for the fermions. 

At large A^ , the integral in (|A.1| ) could be evaluated by the saddle point method. 
The equations determining the saddle points of /x = im and p are 



c sinh(m) + -sinh(— )) = — (A. 2) 



d . .,rn q 
-smhl — I) = — 
2 ^2" p 

and 

TTl 

p{a + cosh(m) + (icosh(— )) = p (A. 3) 

We would like to see weather there is a solution with p = 0. As the right hand side of 
( [A.2|) becomes large in the limit p 0, we can self consistently approximate cosh(m) 
and sinh(m) as e"* and we get 

m ^ log (A. 4) 

cp^ 

Hence a logarithmic potential for p is once again generated. One can also confirm 
this by putting (^) in (pOI) . 



B. Appendix: Positivity of the coefficient of the quadratic 
term in the effective action 

Let us consider the partition function of YM theory on a compact manifold written 
as an integral over the effective action of p = TrUTrU^^ . 

Z{/3) = J D[/e^'(s<=//('')) (B.l) 
= J rfpe^'(^^//('')-^^^(''» (B.2) 

Where Sm{p) is the contribution from the measure part^^ of path integral and 

%;(p) = a(/5)p2 + ^a„(/3)p" (B.3) 

?i=4 

i.e. a polynomial in p. As /? — cxd we have S'e//(p) — > 0. Contribution from the 
measure part S'm(O) has only one minimum at p = 0. Hence at low temperature the 

-^^See discussions before ( |3.20| ). 
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system will have a saddle point at p = 0. Expanding p around this saddle point as 
p = + ^ we get 

Z{(3) = r d{6p)e-^^-''^^^^^^''^'-"'^^^ (B.4) 



Like any thermal partition function, (|B.2|) or ( p.4|) is a decreasing function of the 



p. Hence a(/3) should also be a decreasing function of /3. Since a(oo) = 0, for any 
finite (3, a{(3) is a positive decreasing function. Hagedorn transition happens when 
(^IPh) = 1; but whether a{(3) will reach 1 or not depends on the model. 
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